Quasiperiodicity, bistability and chaos in the Landau-Lifshitz equation 
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The dynamics of an individual magnetic moment is studied through the Landau-Lifshitz equation 
with a periodic driving in the direction perpendicular to the applied field. For fields lower than 
the anisotropy field and small values of the perturbation amplitude we have observed the magnetic 
moment bistability. At intermediate values we have found quasiperiodic bands alternating with 
periodic motion. At even larger values a chaotic regime is found. When the applied field is larger than 
the anisotropy one, the behavior is periodic with quasiperiodic regions. Those appear periodically 
in the amplitude of the oscillating field. Also, even for low values of the driving force, the moment 
is not parallel to the applied field. 
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Traditionally, the study of the dynamics governed by 
the Landau-Lifshitz equation is related to the ferromag- 
netic resonance problemsEJ. Recently, the spin dynamics 
has become also important in other physical phenomena 
relevant to technological .applications, such as, e.g., mag- 
netic recording processes™ due to a continuous increase 
of the magnetic 1 recording density together with the writ- 
ing frequencvtJu . The dynamical micromagnetic calcu- 
lations have provided a useful tool in studying such ira- 
portant media characteristic as dynamical coercitivityEl. 
In spite of the fact that the Landau-Lifshitz equation 
is widely used in micromagnetic calculations, up to our 
knowledge, no systematic study of its dynamics exists in 
the literature. Let us recall here that this equation is 
nonlinear, and in some regime one may expect a highly 
complicated dynamics similar to one arising for an exter- 
nally driven pendulum. As an example, we can mention 
the nonlinear stochastic resonance behavior of an indi- 
vidual magnetic momentO. The purpose of this paper 
is to present a systematic study of nonlinear dynamics 
governed by the Landau-Lifshitz equation, including its 
bifurcation diagram and stability properties. This would 
be very useful when analyzing results obtained from large 
simulations of coupled Landau-Lifshitz (LL) equations. 
There, and for some values of the coupling parameters, 
the individual characteristics of each magnetic moment 
may play an important role in the collective behavior. 
Eventually, when enough of these moments are coupled, 
a description in terms of magnons would be possjJpleEH 

Also the subject of chaos in magnetic materialsQ is not 
new. It has been studied in YIG, both experimental and 
theoretically, through spin wave descriptionstm. In dif- 
ferent driving regime, the chaotic behavior can arise in 
magnetostrictijwe wires and ribbons due to the magnetoe- 
lastic couplingt2l. ._. 

In the original form (Landau-Lifshitz-GilbertE 2 ]) the 
equation may be written as: 

M t = -g (M x H cff ) + (M x M t ) , (1) 
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-M t = — (M x H cff ) — —j- (M x (M x H cff )) . 
9 M 

(2) 

Where g is the local giromagnetic factor, r\ is the damp- 
ing, Mq the saturation magnetization, M is the tridimen- 
sional local continuous magnetization, whose module is 
conserved (M • M = Mo), and H e ff the effective field: 
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(3n (n-M) + aV 2 M + H d . 



(3) 



Here H cxt is the external magnetic field, is the 
anisotropy coefficient, n is the unitary vector pointing 
in the anisotropy direction, a the coefficient of exchange, 
and, Hd the demagnetizing field that takes into account 
the dipolar interactions between local moments in the 
system. 

Here we will deal with only one local magnetic moment, 
so that we will not consider exchange and dipolar inter- 
actions. This description could be relevant tOpdynamics 
of magnetic microwires with a high anisotropj£3. 

The external field will be decomposed in two parts: 
H„, constant and parallel to the anisotropy direction, 
and h = ho ■ sin(wi), perpendicular to n and oscillating 
in time with frequency u. 

For practical purposes we rewrite the equation in adi- 
mcnsional form and expanded notation: 

nm x = - (m y (h z + m z ) - m z h y ) - r\ (m x m z (h z + 

m z ) + m x m y h y + (ml - l)h x ) , (4a) 
Kniy = (m x (h z + m z ) - m z h x ) — tj (m y m z (h z + 

m z ) + m x m y h x + (m y - l)h y ) , (4b) 
nm z = - [m x h y - m y h x ) - r\ (m x m z h x + 
m y m z h y + (m 2 - l)(h z + m z )) , 



or in the more practical form (Landau-Lifshit 



zO) 



(4c) 

where m = M/M , h = H oxt /(/3M ), k = 1 + rj 2 , and 
the dot represents the derivative with respect to an adi- 
mensional time r = g(3Mot. We take the z axis in the 
direction of n and the x axis in the direction of ho , which 
makes h y = 0. 
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These equations with appropriate initial conditions 
have been solved by using a fourth order Runge-Kutta 
scheme. We have observed that the condition that 
|m| = 1 is fulfilled with a precision of more than eight 
orders of magnitude (1CP 8 in 1), for even more than 10 9 
integration steps. 

When looking for chaotic behavior in driven systems, 
the usual quantities to vary are the frequency and the 
amplitude of the perturbation. Sweeping in frequency 
we have found that the interesting ones are those close 
to the resonance frequency, u) r = 1, in our units, and 
of the order of several GHz in real units. Thus possible 
perturbing signals are radio- wave sources. For simplicity 
in what follows will put the perturbing frequency fixed 
to the resonance, and sweep in the amplitude, h = (h^, + 

K) 1/2 

A bifurcation diagram is shown in figure Wfor h z = 0.1 
and r\ = 0.05. Here only the components nig and are 
drawn, though in the numerical simulations we followed 
the equations 4a-c. The diagram shows the values of the 
components of m, at time intervals multiple of T = 2it/uj 
(Poincare sections) , with the value h shown in the x axis 
of the figure. There, various kinds of behaviors can be 
distinguished: periodic, quasiperiodic and chaotic mo- 
tion. When there is only one point for a given value of h, 
it represents a periodic motion with period T; and when 
there is a continuum of points the behavior is quasiperi- 
odic or chaotic. We will now try to describe the principal 
features of different critical points shown in figure |l|. The 
changes in the diagram found when the dissipation (77) is 
changed, are mainly quantitative (changes the value of h 
at which a given critical behavior is found). Also, the dia- 
gram shown has been produced for the magnetic moment 
pointing in the direction of the external field at t = 0. 

For small values of h (h ~ h z ) a discontinuity that 
corresponds to a folding bifurcation is found. This effect, 
consisting of two independent limit circles, may be the ex- 
perimental source for the observation of hysteresis when 
changing the amplitude of the perturbation (h). When 
increasing h the jump in 9 is the one shown in figure [l], 
but if h is decreased, the jump to smaller 9 would occur 
at a lower value of h. Experimentally, this phenomenon 
could manifest itself in a bistable behavior of a magnetic 
microwire near the resonance frequency. 

Two bifurcations, identified as torus have been ob- 
served at h = 0.60 and h = 0.75, leading to two regions 
of complex (but mainly quasiperiodic) behavior: h from 
0.60 to 0.66, and from 0.75 to 0.87. The first one in 
the upper hemisphere, and the second one in the lower. 
In the torus bifurcations the stable limit circles become 
unstable and give rise to quasiperiodic motion on the sur- 
face of a torus. In figure || is shown the phase portrait, 
in coordinates x, y, and z, of the quasiperiodic attractor 
at h = 0.64. There, the Poincare section changes from 
just one single point to a closed connected curve. 

The region from 0.90 to 1.00 is a mixture of periodic 
and quasiperiodic behavior, and even chaotic motion. 
The chaos in this region is characterized by a chaotic 



attractor at h ~ 0.9787, which develops via a global bi- 
furcation of the type of chaotic transients. This means 
that the system will evolve in the chaotic attractor for 
some time, and then, feeling the periodic or quasiperi- 
odic stable orbits, will leave it. This is illustrated in 
figure |^, where the time evolution of the Lyapunov expo- 
nents is shown. Initially both exponents converge, one to 
a positive value, and the other to negative, as signature 
of chaos. But at a given time the positive exponent initi- 
ates a decrease towards 0, or even negative. The Poincare 
sections for the initial and final time steps are also shown. 
Initially the trajectories follow a chaotic map, but after 
some time they eventually fall in a period-seven orbit. 
The time spent in the chaotic behavior becomes larger 
as the chaotic attractor is approached. ■— ■ 

Next, there is a wide region of period doublinglil, with 
some higher period stripes. Finally, from 2.20 to 2.60, 
clear chaotic regions (see the Lyapunov exponents in fig- 
ure [l]) alternating with quasiperiodicity are observed. In 
figure ^ we show the Poincare section of the chaotic at- 
tractor corresponding to h=2.50. In this case the route 
to chaos is also that of chaotic transients. 

If the initial state for the magnetization is in the direc- 
tion opposite to the external bias field, h z , then, basically, 
the picture presented above holds. Nevertheless the fold- 
ing disappears, and the stable period-one orbit evolves in 
the lower hemisphere. The same change of hemisphere 
happens for the period doubling region. The two (upper 
and lower) torus bifurcations are also preserved. 

When the external applied bias, h z , is larger than the 
anisotropy field (larger than 1 in our units) the behavior 
is slightly different. The bifurcation diagram is shown 
in figure [|. The same kind of structures repeats at 
h = 3.25 x n, where n is an integer. A magnification 
of those structures is shown in the inset, where it is seen 
that they consist of a torus bifurcation and several peri- 
odic windows. The magnetization does not remain at a 
value close to the saturation, but wanders over the whole 
sphere. 

In conclusion, we have demonstrated that the dynam- 
ics governed by a driven Landau-Lifshitz equation in cer- 
tain parameter region can be very complicated. We have 
observed manifestations of this behavior: the quasiperi- 
odicity alternating with periodic motion, the transient 
chaotic behavior and finally a chaotic attractor. The ex- 
perimental observation of these predictions is limited to 
the materials which could have a well defined single fer- 
romagnetic resonance frequency, i.e. may be described 
by a single magnetic moment. Nevertheless, some fea- 
tures of this complicated behavior could be persistent in 
a coupled system til The estimation of the parameters 
show that even in the case of a well-defined single mag- 
netic moment behavior, the experimental observation of 
the chaotic behavior would require a powerful source of 
radio-waves. Nevertheless, a folding bifurcation (bistable 
behavior) occurs for parameter values accessible by nor- 
mal radio- wave antennas. 
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FIG. 1. Bifurcation diagram for the 9 and (f> components 
of the magnetization and Lyapunov exponents when h z —0.1 
and n = 0.05. 
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FIG. 2. Quasiperiodic attractor at h — 0.64 in the diagram 
of figure 111 
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FIG. 3. Time evolution for the Lyapunov exponents when 
h = 0.9785 in the diagram of figure 111 The left inset is the 
Poincare section taken from t = to t = 1.8 x 10 5 . The right 
one is taking i = 4.5 x 10 5 to t = 6.3x 10 5 . 



FIG. 5. The same bifurcation as in fig. |l| for h z = 2. Insets 
show the magnification of the region of h between 6 and 7. 



FIG. 4. Chaotic attractor at h = 2.5 in the diagram of 
figure |l]. 
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